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Abstract 

Continuous families of solitons in generalized nonlinear Schodinger 
equations with non-PT-symmetric complex potentials are studied ana¬ 
lytically. Under a weak assumption, it is shown that stationary equations 
for solitons admit a constant of motion if and only if the complex potential 
is of a special form g {x) + ig'{x), where g{x) is an arbitrary real function. 
Using this constant of motion, the second-order complex soliton equation 
is reduced to a new second-order real equation for the amplitude of the 
soliton. From this real soliton equation, a novel perturbation technique 
is employed to show that continuous families of solitons always bifurcate 
out from linear discrete modes in these non-PT-symmetric complex po¬ 
tentials. All analytical results are corroborated by numerical examples. 


1 Introduction 

Nonlinear wave systems fall into two major categories: conservative and dissi¬ 
pative. Conservative systems are energy-conserving, and their solitary waves 
(solitons) exist as continuous families with continuous ranges of energy values. 
A typical example is the nonlinear Schrddinger (NLS) equation. Dissipative sys¬ 
tems contain gain and loss, and their solitons are generally isolated with certain 
discrete energy values. A typical example of this type is the Ginzburg-Landau 
equation. A recent discovery is that, in dissipative but parity-time {'P'T) sym¬ 
metric systems, solitons can still exist as continuous families with continuous 
energy values [il[lJ|lll[3ll[71IHl[SJ[ini[IIl[Il[Il[Il[I3[Ill[IIl[IHl[Il. An 
example in this category is the NLS equation with a complex but V'T -symmetric 
potential. These soliton families are allowed since the VT symmetry assures that 
the gain and loss of the soliton is perfectly balanced at arbitrary energy levels. 

In dissipative and non-UT-symmetric systems, the expectation is that any 
solitons will be isolated with discrete energy values, as seen in typical dissipative 
systems |20| . However, exceptions were reported numerically in for the 
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NLS equation with a non-PT-symmetric complex potential of special form, 
where families of solitons with continuous energy values can bifurcate out from 
the linear modes of the potential. This finding is very surprising in view of the 
lack of T^Tsymmetry here. For these special potentials, a constant of motion 
was discovered in [52] for the stationary soliton equation. Using this constant of 
motion, soliton families in these special potentials were explained by a numerical 
shooting argument [22] . 

In this article, we analytically investigate solitons of the NLS equation with 
non-PX-symmetric complex potentials. We focus on three main questions; (1) 
What types of non-T^T-symmetric complex potentials admit soliton families? 
(2) How can one analytically explain and calculate soliton families bifurcating 
from linear modes in such potentials? (3) Do these soliton families exist under 
other nonlinearities? 

Regarding the first question, we recognize that in the absence of V'T symmetry, 
the existence of a constant of motion in the stationary soliton equation plays a 
crucial role in the existence of soliton families. Assuming this constant of mo¬ 
tion for complex potentials is a continuous deformation of one that exists in the 
NLS equation without a potential, we show that the only complex potentials 
which admit a constant of motion are those of the form reported in [211 [22], i.e., 
V{x) = g^{x) + ig'{x), where g{x) is an arbitrary real function. This strongly 
suggests that potentials of the above form are the only one-dimensional non- 
T^T-symmetric complex potentials that admit soliton families. 

On the second question, through use of the constant of motion, we reduce 
the second-order complex soliton equation to a new second-order real equation 
for the square of amplitude of the soliton, which is then solved perturbatively 
for a continuous range of g, values. This way, the existence of soliton families 
bifurcating from linear modes in non-T^T-symmetric potentials is analytically 
explained and explicitly calculated. Interestingly, this perturbation calculation 
of solitons differs significantly from the method used for real and V'T -symmetric 
potentials because the linearization operator of the new real equation has a 
distinctly different kernel structure. 

Regarding the third question, we show that these soliton families still exist 
under a more general class of nonlinearities. Furthermore, the choice of nonlin¬ 
earity within this class has no effect on the existence of a constant of motion. 

These analytical results are compared with numerical examples, and good 
agreement between them is illustrated. 

2 Preliminaries 

The mathematical model we consider, in most parts of this article, is the NLS 
equation with a complex potential, 

*^'t + «'xx + R(a;)«' + cr|Tp^- = 0, (1) 

where cr = ±1 is the sign of cubic nonlinearity. This model describes paraxial 
nonlinear light propagation in a waveguide with gain and loss [2] [23] , as well 
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as Bose-Einstein condensates with atoms injected into one part of the potential 
and removed from another part of the potential [24l [25]. Most of the earlier 
work focused on the case where the complex potential V{x) is PT-symmetric, 
i.e., V*{x) = V{—x)^ with the superscript representing complex conjugation 
nn Ilia niii in]- In this article, we consider the case where V{x) is not 
T^T-symmetric, i.e., 

V*{x) V{-x). (2) 

Soliton solutions of equation m take the form 

(3) 

where tjj{x) is a localized function solving the stationary equation 

il^xx-^J'tp + V{x)tj; + a\'^|J\'^1|J = 0, (4) 

and ^ is a real propagation constant. In PT-symmetric potentials, solitons 
exist as continuous families parameterized by /r 01 uniiHiin]. But in non- 
T^T-symmetric potentials, soliton families are generically forbidden [20]. Sur¬ 
prisingly, it was reported recently through numerical examples that in complex 
potentials of the special form 

V{x) = g'^ix)+ig'{x), (5) 

where g{x) is an arbitrary real function, soliton families can still bifurcate out 
from linear modes even when V{x) is non-T^T-symmetric (i.e., when g(x) is not 
even ) [2l[22]. This result is very unintuitive. Indeed, if one performs a reg¬ 
ular perturbation calculation of soliton families bifurcating from linear modes 
in a general complex potential, it will be seen that infinitely many nontrivial 
conditions would have to be satisfied simultaneously, which makes such bifur¬ 
cation almost impossible [20] • However, for the special complex potential (0), 
all those conditions are met, which is miraculous. Obviously this phenomenon 
needs better understanding. A step in this direction was made in |22j . where 
through the discovery of a constant of motion for the soliton equation (0) under 
the potential (0), soliton families in Eq. (0) were explained through a numerical 
shooting argument. 

Many important questions are currently open regarding soliton families in 
non-PT-symmetric complex potentials. For instance, what other non-PT- 
symmetric complex potentials admit soliton families? How can one analytically 
explain and explicitly calculate soliton families bifurcating from linear modes in 
non-PT-symmetric potentials? Do these soliton families also exist under other 
nonlinearities? These questions will be investigated in the remainder of this 
article. 


3 Constant of motion 

A quantity J{x, tp) is called a constant of motion in the stationary equation 0) 
if dJ/dx = 0. The existence of a constant of motion proves to be important for 
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the existence of soliton families (see [52] and later text). Thus, in this section, 
we study what complex potentials V(x) admit a constant of motion. In this 
study, solutions, tpix), to the stationary equation (|T|) are allowed to be any 
solutions, not necessarily solitons. That is, ip(^) is allowed to be non-local. 

First we split the complex potential V(x) into real and imaginary parts, 

V(x) = vi(x) + iv 2 (x), (6) 

where vi(x), V 2 (x) are real functions. We also express the complex function 'i/’(x) 
in polar forms, 

V'(x) (7) 

where r(x), 0{x) are real amplitude and phase functions. Substituting these 
expressions into the soliton equation (|3|), we get 

rxx — fJ-r + vir + ar^ — 9^r = 0 , ( 8 ) 

{r^d)x = -V2r^- ( 9 ) 

In the absence of the potential (z;i = U 2 = 0), it is easy to verify that this 
system admits two constants of motion 

Ji = r^O, ( 10 ) 

and 

J, = rl-fir^ + ^r^ + r^e\ (II) 

where dJi/dx = dJ 2 /dx = 0. Since the system is third order, these are the only 
constants of motion the system can allow. These two constants of motion are 
associated with the flux terms of the power and momentum conservation laws of 
the potential-free NTS equation, but this fact is not important to our analysis. 

In the presence of the potential, it is reasonable to assume that the cor¬ 
responding constant of motion J is a continuous deformation of those in the 
potential-free case. In other words, J approaches constants of motion of the 
potential-free equation when vi,V 2 approach zero. Notice that Ji and J 2 have 
different ranks [26] . Thus under the limit vi,V 2 —t 0, J can only approach 
one of (Ji, J 2 ), not their linear combination. Our strategy then is to calcu¬ 
late dJk/dx (fc = 1, 2) in the presence of the potential and derive conditions on 
{vi,V2) SO that dJk/dx is a total derivative of x, i.e., a constant of motion is 
admitted. 

Now we calculate dJk/dx in the presence of a potential. First we consider 
dJi/dx. Eq. (O clearly shows that, in order for dJi/dx to be a total derivative, 
we must have V 2 = 0, i.e., the potential V{x) is real. This is not what we want 
since we exclusively consider complex potentials in this paper. Thus, there are 
no constants of motion in Eqs. that approach Ji when the complex 

potential approaches zero. 

Next we consider dJ 2 /dx. Utilizing equations ®-®, we readily find that 

^ =-vi{r'^)x-‘2v2r'^0. ( 12 ) 
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The right side of this equation can be rewritten as 


7 T n 

= Wx + r^vij,+2[r'^e)j, I V2dx, (13) 

where 

W = —vir^ — 2r‘^6 J V2dx. 

Then utilizing equation ([9]), the above equation becomes 

^ =Wa:+r'^ (^Via; - 2v2 J 

In order for the right side of the above equation to be a total derivative, the 
necessary and sufficient condition is 



Vix = 2 v 2 j V 2 dx. 
This condition can be rewritten as 

V2dx 


Vlx = 

thus ^ 

vi = V2dx^ + C, 

where C is an arbitrary constant. Finally, denoting 

g= V2dx, 


(15) 


(16) 


(17) 


(18) 


the potential V(x) which admits a constant of motion then is of the form 

V(x) = g^(x)+ig'(x)+ C. (19) 

Obviously, the constant C in this potential can be eliminated from Eq. m 
through a simple gauge transformation. The remaining potential is then of the 
form Thus we conclude that if the constant of motion for the stationary 
equation (|4]) with a complex potential is a continuous deformation of J 2 without 
the potential, then this constant of motion exists if and only if the complex 
potential V(x) is of the special form ([S]), and the corresponding motion constant 
is 

J = J^-W = J 2 + g^r^ + 2gr'^e, ( 20 ) 

or more explicitly, 

J = rl-^,r^ + y^ + r^{e + g)\ ( 21 ) 

where dJ/dx = 0. This constant of motion agrees with that reported in [22] for 
these special potentials ©■ 
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4 Bifurcation of soliton families 


In this section, we analytically calculate the bifurcation of solitons from linear 
modes in Eq. o. with potential of the special form (O, and show that soliton 
families bifurcate out in such non-PT-symmetric systems. 

For the potential ([S]) , when solitons ([3]) are expressed in polar forms (O , the 
equations for r and 9 are seen from Eqs. (0)-® as 

Txx - ^ir + g^r + or^ - 6>^r = 0, (22) 

{r^9)a: = -gxT^, (23) 

and these equations admit a constant of motion (EH. For solitons, this constant 
J can be evaluated at a; = oo as zero, thus 

+r^{6 + gY = 0. (24) 

From this equation, we get 

e = -g± J. (25) 

Inserting it into Eq. (l22|) and after simple algebra, we get 


rx:x - 2^r + -ar + y = ^2^y 9-r “ 


(26) 


or 

(r^)xx - 4^r^ + 3 ot'^ ± 2g\/4/xr^—= 0. (27) 

Denoting R = , we arrive at a single second-order equation for the real am¬ 

plitude function R as 

Rxx - 4/ii? -f 3cri?2 ± 2gy/AgR'^ - 2aR^ - Rl = 0, (28) 


which can also be rewritten as 


_ 2aR^ - Rlj = ±2gR^. 


(29) 


4.1 Perturbation calculations 

The sign in Eq. (|28ll needs to be chosen appropriately according to the function 
g{x). Indeed, if g{x) switches to —g{x), this sign should switch as well. Without 
loss of generality, we take the plus sign in Eq. (l28|) . 

Rxx - 4/ri? -f 3cri?2 -b 2g^/AgR'^ - 2aR^ - R^ = 0. (30) 

Note that sometimes the same solution R{x) can lead to mixed signs in Eq. (I28|) 
on different x-intervals. This could occur if AgR^ — 2aR^ — R^ is zero somewhere 
on the a:-axis, since the square root is a possible mechanism for inducing a sign 
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change so that the square-rooted quantity remains smooth. We do not consider 
such mixed cases here. This exclusion will be assured by Assumption 1 in 
Sec. 14.31 

For a localized function g{x)^ it is easy to see that the large-a: asymptotics 
of the soliton solution in Eq. (1301) are, to leading order, 

i?(x) —)► X —>• ±oo, (31) 

where a± are positive constants. 

First we consider linear modes in Eq. O, which satisfy the equation 

(t^xx - Vo<(> + 25-v/4/ro<()2 - 4>l = 0. (32) 


This equation can be rewritten as 

(\/4Aio<(>^ - = ‘^9(t>x- (33) 

Since g and (j) are localized functions, we see that 

/ X noc 

2g(j)^d^ = - 2g(j)^d^, (34) 

-oo J X 


and 


= 0 . 


(35) 


Equation (1321) is scaling-invariant and thus an eigenvalue problem, but it is 
nonlinear in both the eigenvalue go and eigenfunction (j). Thus, this is a different 
type of eigenvalue problem. Solving this new eigenvalue problem is equivalent 
to solving for discrete real eigenmodes in the original eigenvalue problem from 
Eq. dH), i.e., 

i’xx + V {x)iIj = goip, (36) 

and the eigenfunction correspondence is 4> = |'0p. Previous results in [21] have 
shown that for the underlying special potential ([5]) , the linear eigenvalue problem 
(l36l) admits discrete real eigenvalues for a large class of functions g{x). The new 
eigenvalue problem (1321) makes the existence of such real eigenvalues more clear 
since all quantities in that equation are real. 

From such eigenmodes, families of solitons can bifurcate out under varia¬ 
tion of g. We will analytically prove this by explicitly calculating this soliton 
bifurcation from a linear mode {go, <p) using perturbation methods. 

The perturbation expansion is 


R — ^{Ro T • )i (3f) 

g = go + €, (38) 


where e > 0 is a small parameter. Here we have assumed the bifurcation occurs 
to the right side of go- As we will see later, this assumption dictates the sign of 
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nonlinearity a. If the bifurcation occurs to the left side of /xq, then only trivial 
modifications to our analysis are needed, and the bifurcation will occur for the 
opposite sign of nonlinearity. 

Inserting the above expansion into Eq. o, at order e, we hnd 

Ro = co(j), (39) 


where cq is a positive constant to be determined. 
At order e^, we get 

LRi = F, 


where 


Pi = 


L — dxx + Pidx + P2, 

‘^94’x . I ‘^94’ 




P 2 = 4//0 


A/4/ro(/)2 - (/)2 


- I 


and 


fl=4:(j)\ 1 - 


94^ 


y'4:fio4>'^ - 4>l 


F = co(/i - cocr/ 2 ), 


/2 = 3 - 


294) 


V^4p(o(/)2 - (pl 


(40) 

(41) 

(42) 

(43) 

(44) 


Now it is time to analyze the properties of homogeneous solutions and ad¬ 
joint homogeneous solutions of the operator L and the solvability condition of 
Eq. dlOl. 


4.2 Kernels of linearization operators L and 

First, it is easy to verify that ^ is a homogeneous solution of L, i.e., 

L(j) = 0. (45) 

Let us suppose the other homogeneous solution of L is (j) 2 , then according to 
Abel’s formula, the Wronskian of (</), ^ 2 ) is 

W{cj),^2) = Woe-fP^‘^^, (46) 

or 

Will), (j)2) = lEox/W^^ (47) 

in view of Eqs. (l33l) and (1^ . Here Wq is a constant. Utilizing Eq. (l34l) . the 
above Wronskian can be rewritten as 

/ X noc 

294)^d^ = -Wo / 2g^^dC (48) 

-00 J X 

From this formula we see that, if (/(x) is a localized function, then the decay 
rate of this Wronskian at large |x| is faster than that of (j), and thus ^2 is also 
a localized function. 
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Using these homogeneous solutions of L, we can build homogeneous solutions 
of the adjoint operator L^, where 

= da:a: “ + P2- (49) 

Lemma 1 The two homogeneous solutions (0^, (j)^) of adjoint operator 


are 


02 




yj 

fU(0,02)’ ~ 

Proof: We first turn the second-order homogeneous equation of operator L 
into a system of first-order equations, 


X' = 


0 1 
-P2 -Pi 


X, 


(51) 


where the prime stands for derivative to x. The fundamental matrix solution 
to this system is 

0 02 
0a: 02x 


X = 

The adjoint system of Eq. m is 

r' = - 


(52) 


0 -P2 

1 -pi 


Y. 


(53) 


Notice that if U = [yi,y 2 ]'^, where the superscript ‘T’ represents vector or 
matrix transpose, then it is easy to verify that 


L^y2 = 0 , 


(54) 


i.e., the second component of vector solution Y is in the kernel of the adjoint 
operator L^. 

It is well known that the fundamental matrix solution to the adjoint vector 
system (l5^ is (X“^) . This can be proved by calculating (XX“^)', where 
upon utilizing Eq. (EH, a homogeneous differential equation for X ^ would 
be obtained. Taking the transpose of this equation would reveal that 
satisfies the adjoint equation (l53)) - Notice that 

^ 

^ ’ W^(0,02) 

Since the second-row functions in this matrix are in the kernel of the adjoint 
operator the functions 0"^ and 0^ defined in Eq. (l50l) are then homogeneous 
solutions of the adjoint operator L^. □ 

In view of Lemma 1, if g{x) is a localized function, then both adjoint homo¬ 
geneous solutions 0"^,0^ are unbounded, because the decay rates of 0 and 02 
at large |a:| are slower than that of the Wronskian IU(0, 02). 

The fact that L has only localized solutions and has only unbounded solu¬ 
tions in their kernels makes the solvability condition for the first-order equation 
((4(1)) novel, as we will delineate below. 
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4.3 Solvability conditions for certain potentials 

In this subsection, we show how to impose the solvability condition on Eq. (HOI) 
under the following assumptions. 

Assumption 1 For the linear eigenmode ~ 4>x strictly 

positive for all x] 

Assumption 2 The function g{x) decays exponentially at large x as 

gix) ^ , a; —>■ ±oo, (56) 

where b± and /? > 0 are constants. 

Assumption 3 For these potentials, > f3. 

Assumption 1 assures that the linear operators L, are nonsingular. In 
addition, there will be no sign change on the cc-interval in Eq. (l28ll . This 
assumption will be made throughout the text. 

Assumptions 2 and 3 are introduced in order to make our analysis more 
explicit. If the function g{x) does not satisfy these assumptions, an alternative 
analysis will be outlined in the next subsection. 

Remark 1 In some sense Assumptions 2 and 3 represent the most common 
case, since for practical purposes the exact decay rate at large x should have 
minimal effect on the dynamics of the system. Hence modifying the small tails 
of the potential to have a suitably exponentially decaying rate should not make 
a meaningful difference. However, we will still show in the next subsection that 
these perturbation calculations may be performed for general potentials, and 
then verify all results numerically. 

At large |x|, the asymptotics of the eigenfunction 4>{x) can be readily seen 
from Eq. (1321) as 

(/)(x) —>■ a; —>■ ±oo, (57) 

where ')± are constants. Then under Assumption 2, it is easy to see from Eqs. 
((Ml) that the large-x asymptotics of _ 02 jg 



\/4/ro(/>^ — 0 ^ —t s±e x —>■ ±oo. 

(58) 

where 

25±7±i/4^ 



/3 + v^ ■ 


Thus 

Pi(x)-?■ ±(/3 + V?^), x->-±oo. 

(59) 


P 2 {x) —>■ /3V4/ro, X —>■ ±oo. 

(60) 

hence the asymptotics of operators L and are 



da:x±{l3+ \/4^)9x + /3^/4^, x ±oo. 

(61) 

and 

t dxx T (/3 + \/ ^go)dx + /3\/4/ro, x —>■ ±oo. 

(62) 
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From these asymptotics, it is seen more explicitly that all homogeneous solutions 
of L are localized (as or their linear combinations), and all 

homogeneous solutions of are unbounded (as , or their linear 

combinations). 

Regarding the second homogeneous solution 4>2{x), in view of the asymp¬ 
totics (IS3 of the first homogeneous solution (^(x), without loss of generality we 
can set the large negative-x asymptotics of 4 > 2 (x) as 

02(x) —>■ X —>■ — oo. (63) 

Then its large positive-x asymptotics is 

4 ' 2 { x ) —>• -I- X —>■ -boo, (64) 

where ki, H 2 are constants. Substituting these asymptotics into the Wronskian 
function and using the Wronskian formula (H51) . the value of ki can be deter¬ 
mined. But this K\ value is not needed in our analysis. 

From Lemma 1 and Eq. (SB, we rewrite the adjoint homogeneous solutions 
(j)'^ and (j )2 equivalently as 


\/4/ro<?5>2 - (j)l 


^2 




(65) 


Then using the asymptotics (1571) . (1551) . (1551) and (IMl) . we find that the large-x 
asymptotics of (j)^ and are 




X —>• —oo. 



^ , X — >• -|-00, 

(66) 

and 

(t>2{x) -t 

X — >• ±oo. 

(67) 

where 

^ _ /3 + 

26±7±i/47to’ 

J± = d±j±. 

(68) 


Notice that both adjoint solutions are unbounded and grow exponentially at 
large x. 

The asymptotics of functions fi and /2 in the first-order equation (|40)) can 
be similarly obtained as 



/i(a^)^<Z±e-^l"l, 

X —>■ ±oo. 

(69) 

and 

/ 2 (x) ^ 

X —>■ ±oo. 

(70) 

where 


II 

I 

^ I 

-.J 

H-io 

(71) 
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Now we consider the solvability condition of the first-order equation (HUl) . 
We see from Eqs. (1^ . (l40l) . (1^ and (TfOl) that the large-a; asymptotics of Ri{x) 
must be 


Ri{x) -b x ±oo, (72) 

where (x) are certain linear functions of x (these linear functions come about 
when one expands the tail function of (ISTI) into a perturbation series 

around /i = /tq), and are constants. Note that the tails in the 

above equation are induced by the nonlinearity-related forcing term /2 and are 
admissible. They do not contradict the leading-order R{x) asymptotics (1311) 
since they are of higher order. Enforcement of this tail behavior for Ri{x) will 
yield the solvability condition which determines the cq value. 

We begin by taking the inner product of Eq. (l40l) with to get 

= ((/)^,Co(/i - Cocr/ 2 )), (73) 

where the inner product is defined as 



Performing integration by parts, the left side of this equation becomes 
{cl>tLRi) = (T^0^,i?i)+ 

= i^^Ri.-cl>iRi+Pi<i>^Ri)CZ- ( 75 ) 

In view of the asymptotics of (j)^, Ri and pi in Eqs. (EZl) and dZlD, as well 
as Assumption 3, we see that the right side of the above equation is zero, hence 
we obtain a solvability condition from Eq. dzi as 

- C 0 CT/ 2 ) = 0. (76) 


This solvability condition is the analog of Fredholm Alternatives condition, and 
it quickly yields the formula for cq as 


ct(<(>2 ,/2) ’ 


(77) 


Notice that /i and /2 decay at large a; as e or faster [see Eqs. (IMll - dTOl) ]. 

and 4 > 2 {^) grows at large x as [see (EzH- Thus under Assumption 3, both 
integrals in the inner products of the above equation converge, and hence cq is 
well defined. 

Eq. (1771) is a necessary condition for the existence of the first-order solution 
Ri{x) with suitable asymptotics (1771) . Since cq must be positive, Eq. (1771) then 
shows that, in order for the soliton bifurcation to occur to the right side of po 
[see Eq. (l38l) ]. the sign of nonlinearity tr must be chosen as the sign of the ratio 
{(l>2 Jl)/{(t^2 J 2 ) ■ 
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The above solvability condition (EZl) turns out to be also sufficient for the 
existence of solution Ri{x) with suitable asymptotics (m. To show this, we 
notice that the general solution to the first-order equation dini) can be derived 
by variation of parameters as 




h{x) [ 
Jo 


(t}2i0F{0dC + Cicl){x) + C2(l)2ix) 


. 

where ci and C 2 are real constants. Using the asymptotics detailed earlier in this 
section and under Assumption 3, we find that at large x, (j)^{x)F{x) approaches 
a constant, and (^^(x)F(x) decays exponentially. Thus 


Jo 


X —>■ zbcxi, 


(79) 


where P^{x) are linear functions of x, and 


px ^±00 

/ / <l> 2 {£)F{ 0 d^. x^ioo. 

JO JO 


( 80 ) 


In view of the large-x asymptotics of ^(x) and 4>2(x), in order for i?i(x) in (I78|) 
to exhibit the suitable asymptotics (1721) . the necessary and sufficient conditions 
are 


^±00 

/ </>^(C)i"(e)d^ + C2 =0, 

^0 


( 81 ) 


which leads to the equation 


f*+oo 


4>2{x)F{x)dx = 0 . 


(82) 


Substituting the expression (I43p for Fix) into this equation, we then obtain the 
Co formula (EZl). Hence this cq formula is a necessary and sufficient condition 
for the existence of solution i?i(x) with suitable asymptotics (1721) . 

In the i?i formula (l78l) . while co is given by formula (1771) and C 2 given by 
equation (1811) . ci is still a free parameter. This ci parameter will be fixed by 
requiring the second-order solution i? 2 (x) to have suitable large-x asymptotics 
[similar to (17^ but with linear functions (x) replaced by quadratic functions 
H^(x)]. This calculation of ci is in the same spirit of the cq calculation, thus 
details will not be pursued in this article. 


4.4 Extension to general potentials 

In the event that Assumptions 2 and 3 of the previous subsection do not hold, 
i.e., the decay rate of the potential is not simply exponential, or the exponential 
decay rate is too fast, i.e. /3 > then the simple cq formula (l77)) in the 

previous subsection will be invalid. For instance, when /3 > i/d/xo, the integral 
in the numerator of (1771) would be divergent in view of the asymptotics of its 
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integrand. The quantity on the right side of Eq. (75) would not vanish either. 
Thus the solvability condition for these more general potentials needs a different 
treatment. 

In our new treatment, we consider the Ri{x) solution (1781) and demand that 
its tail asymptotics match (1721) . In particular, this entails choosing cq such 
that the terms of which decay like g{x) for exponential potentials) are 

eliminated. 

Suppose the tail asymptotics of the second homogenous solution 4’2{x) has 
the form 


<('2 -t 


r (x), X —>• —oo, 

+ T+(x), X —>■+00, 


(83) 


where K 2 is a certain real constant, and t^(x) are the other decaying tail func¬ 
tions. This (j )2 asymptotics is the counterpart of Eqs. (I63ll - (l64l) in the previous 
subsection. 

We substitute the F formula (H51) into (1751) . Then this Ri solution can be 
rewritten as 


i?i(x) = Co [i?ii(x) - Cocri?i 2 (a:)] -b Cicj){x) + C 2 (l) 2 {x), (84) 

where i?ii(x) and i?i 2 (x) are particular solutions of equations 

Ti7ii(x)=/i, LRi 2 {x) = f 2 . (85) 

For definiteness, we impose zero initial conditions on Rn and i?i 2 at x = 0, i.e., 

i?ife(0) = i?;fe(0) =0, fc = l,2. (86) 

Notice that both particular solutions i?ii and R 12 approach zero at large x, 
since the forcing terms fi and /2 approach zero, and the homogeneous solutions 
are all localized. 

The tails of the particular solutions i?ii and i?i 2 each have terms that decay 
exponentially and a term which decays like t^(x), due to the exponentially 
decaying forcing terms and exponential tails inside the homogeneous solutions. 


Specifically, 




i?ii(x) 

- i 

f P^^{x)e''^^^ + Xi T (x), X —>■ — 00 , 

( P]^(x)e“'/^“-b X)''t+(x), X —>■-boo, 

(87a) 

i?i 2 (x) 

- ] 

f P-^ + Cf -H x"T" (x), 

( p+ x+T+{x), 

X ^ - 00 , 

X —>• -boo. 

Here are 

linear functions and P^, are constants. 



Substituting these asymptotics into the i?i(x) formula (I84|) and comparing 
its tails with Eq. (ESI), we see that the coefficients on t^(x) must be zero as 
X —>■ ± 00 . This leads to the following system of equations 


CoXi - clcrx2 +C2 = 0, coxt - clcrxt +C2=0. (88) 
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From these, we obtain the necessary and sufficient solvability condition as 


Co 


xt -Xi 

^ {xt - X2 ) 


( 89 ) 


For potentials with exponential decay rates, the constants xt^xt can be 
found analytically with a bit of effort. However, in general these constants 
in the tails of Rii{x) and Ri 2 {x) may not be known analytically, since the 
tail behaviors t^{x) of the second homogeneous solution (j) 2 {x) may not be 
analytically available. Regardless, these constants xf and xt can be efficiently 
evaluated numerically, as examples in the next subsection will show. 


4.5 Numerical examples 

Now, we numerically confirm the above analysis with two examples. 

Example 1 For the first example, we choose the complex potential ([S]) with 
an uneven double-hump function 

g{x) = 0.8 [sech(x -I- 2) -|- /isech(x — 2)], (90) 

where h is a positive constant, and a = 1 (focusing nonlinearity). Notice that 
this potential is exponentially decaying, satisfying our Assumption 2 with /? = 1. 

When h = 1.2, the function g{x) and the corresponding complex potential 
V{x) are displayed in Fig. [Ila,b) respectively. Notice that this potential is 
non-PT-symmetric. Eigenvalues of the eigenmode problem (1321) are shown in 
panel (c), where two real eigenvalues are found. The larger of these eigenvalues 
is /io ~ 0.3708, whose eigenfunction cj){x) is plotted in panel (d). For this 
eigenvalue, v^4/ro > 1, thus Assumption 3 is met, and the analysis in Sec. 14.31 
applies. 

From this linear eigenmode, we have verified numerically that a continuous 
family of solitons bifurcates out. The power curve of this soliton family is shown 
in panel (e). Here the power is defined as P = \tjj\‘^dx. The analytical 
prediction for the power slope P'{no) at the bifurcation point can be obtained 
from equations (I571) - (IM1) as 


-Panal(Mo) = Co / (jjdx, (91) 

J —OO 

where cq is given by formula dZTl). For h = 1.2, this analytical power slope is 
found to be approximately 5.8961. The line with this power slope is plotted 
as dashed red line in panel (e), and good agreement with the numerical power 
slope can be seen. In panel (f), the amplitude profile R = IV'P of the soliton at 
the marked point of the power curve (with g, = 0.6) is displayed. 

As parameter h in the g{x) function (1901) varies, the discrete eigenvalue go 
will change [see Fig. Ha)]. When h drops below 0.926, go will fall under 0.25, 
entering the i/4/io < P regime (where Assumption 3 does not hold). In order 
to test our theory for both ^JTgo > P and < P cases, we have plotted in 
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Figure 1: (a) The example function g{x) in Eq. (IM)l with h = 1.2; (b) Complex 
potential V(x) for the g(x) function in (a); (c) Eigenvalues go of the eigenmode 
problem (1321) : (d) Eigenfunction (j) of the largest eigenvalue in (c); (e) Power 
curve of solitons bifurcating from the linear mode in (d) under focusing nonlin¬ 
earity [solid blue: numerical values; dashed red: the line with the theoretical 
slope value (lOTl) ]: (f) Soliton (IV'P) at the marked point of the power curve. 
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Fig. [2Ib) the theoretical predictions for the power slope ^anal(Mo) m Eq. (gD) 
for 0.5 < h < 1.5, which encompasses both cases. The reader is reminded that 
the Co formula is given by Eq. (1771) when > P and by Eq. (IM)l when 

^ /?. In the same figure, numerically obtained power slopes F*^(/ro) for 
each h value are shown as well. It is seen that numerical and analytical slope 
values exactly match each other, confirming the accuracy of our theoretical 
analysis in sections 14.3114.41 



Eigure 2: (a) Linear eigenvalue in Eq. (l32ll for the g[x) function (1^ at 
various h values; the horizontal dashed line is at /tq = 1/4; (b) comparison 
on the slope of the power curve at the bifurcation point between numerical 
values (solid blue) and analytical predictions (red dots). 

Example 2 As for the second example, we consider the potential ([5|) with 

g{x) = (92) 

The resulting potential V{x) is displayed in Eig.[3Ka). The tails of this potential 
decay like a Guassian, which is faster than exponential. Thus the results in 
Sec. 14.41 apply. In this case, analytical expressions for the tail functions t^{x) 
of in Eq. (1831) are not easy to obtain, but numerical approximations can 

be readily computed. Specifically we select the (f )2 (x) function by requiring that 
for X ^ —oo the function decay like a gaussian and take this tail to be t~{x). 
Now for X —>• + 00 , the dominant decay of this tail is exponential, i.e., the tail 
term t^{x) decays faster than in Eq. (l83l) . thus one must first find 

the coefficient K 2 of the exponential tail from large-a; values of (j> 2 {x). Then 
subtracting away this exponential tail from 1 / 2 ( 2 :), the remaining tail is then 
r+(a;). To obtain and X 2 ia Eq. (1871) . we first compute i?ii(a:) and Ili 2 (x) 
from the inhomogeneous equation (HOI) , with F replaced by fi and / 2 , under the 
initial conditions (I86|) . This is done by integrating the inhomogeneous equation 
from X = 0 out to x = ± 00 . By substracting their (slower-decaying) exponential 
tails and comparing the remaining tails with (x) in (f> 2 , X? and x^ can then 
be ascertained. From these numbers, the cq value is calculated from formula 

dMD- 
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Now we compare these analytical predictions against numerical results. Solv¬ 
ing the eigenvalue problem (l32l) . we find three discrete real eigenvalues, the 
largest being /ig ~ 2.6923. From this eigenmode, we have confirmed that a 
soliton family indeed bifurcates out. If the nonlinearity is focusing (tr = 1), the 
power curve of this soliton family is plotted in Fig. Ellb), and the profile of the 
soliton at the marked point of the power curve (with ^ = 3.5) is illustrated in 
panel (c). On the power curve, the line with analytically predicted power slope 
at the bifurcation point from Eqs. (15^ and m is also plotted. It is seen that 
this analytical power slope matches the numerical one very well. 



Figure 3: (a) Complex potential V{x) for the g{x) function in Example 2; (b) 
Power curve of solitons under focusing nonlinearity [solid blue: numerical values; 
dashed red: the line with the theoretical slope (IMT) ]: (c) Soliton (I'f/’P) at the 
marked point of the power curve. 


5 Extension to more general nonlinearities 

In this section, we show that the results in the previous sections can be readily 
extended to a wider class of nonlinearities 

+ + + = 0, (93) 

where G(-) is an arbitrary real function, and V{x) is a complex potential. Soli¬ 
tons (O in this equation satisfy the stationary equation 

tpxx - gip + v{x)tp + GdV'n?/’ = 0. (94) 

Just as in the case of cubic nonlinearity, in the absence of the potential 
\V(x) = 0], this soliton equation admits two constants of motion. Assuming 
that the constant of motion in the presence of the complex potential V (cc) is a 
continuous deformation of those without the potential, we can show by the same 
technique employed in Sec. [3] that the only complex potentials which admit a 
constant of motion are those in the special form of ([5]), and the corresponding 
constant of motion is 


J = rl-gr^+H{r'^)+r^{e + gf, (95) 
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where H{z) = f G{z)dz, and dJ/dx = 0. 

We can also show that for these general nonlinearities, with potentials of the 
form (O, continuous families of solitons still bifurcate out from linear discrete 
eigenmodes. Without loss of generality, we require H{Q) = 0. Then for solitons, 
J = 0. Using this relation, the equation for the complex soliton 'ip{x) is reduced 
to the following second-order equation for the real amplitude variable R{x) = 
\'ip{x)\^: 

R^^-AnR + 2RG{R) + 2H{R)±2gyiiI^^^lRW(R)^^ = 0. (96) 

This equation is the analog of Eq. (1^51) for the cubic NTS equation ([T]). Re¬ 
peating similar analysis as in the earlier text, these soliton bifurcations can be 
explicitly calculated. 

To illustrate these analytical results for general nonlinearities, we consider 
the following example with a saturable nonlinearity. 

Example 3 Consider the NLS equation (IMl) with a saturable nonlinearity 
and complex potential. 


+ + + = (97) 

where the potential V(x) is of the special form ([5]) with g{x) chosen the same as 
in Example 1 [i.e., g{x) is given by Eq. (TO!) ], with h fixed as h = 1.2. Solitons 
in this equation are sought of the form ([3]), where tpi^) solves 

I'd’P 

- g.tp + y{x)ip + ^ = 0. (98) 

Since the potential here is the same as that in Example 1, discrete eigenvalues 
in the linear equation (1321) remain the same as those shown in Fig. [TJc), with 
the larger one being /ig ~ 0.3708. From this eigenmode, we have confirmed that 
a continuous family of solitons bifurcates out, whose power curve is displayed in 
Fig-HDa). At the marked point of the power curve, the corresponding soliton is 
plotted in Fig.lUb). This example verifies that the bifurcation of soliton families 
in complex potentials ([5]) occurs for a wider class of nonlinearities (l93l) . 

6 Summary and discussion 

In this paper, we have analyzed soliton families in NLS-type equations with non- 
T^T-symmetric complex potentials. Under a weak assumption, we have shown 
that stationary forms of these equations admit a constant of motion if and only if 
the complex potential is of the special form ([5|). Using this constant of motion, 
we reduced the second-order complex soliton equation to a new second-order 
real equation for the amplitude of the soliton. From this new soliton equation, 
we showed, by perturbation methods, that continuous families of solitons always 
bifurcate out from linear eigenmodes for this special form of complex potentials. 
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Figure 4: (a) Power curve of solitons bifurcating from a linear mode in Exam¬ 
ple 3; (b) Soliton (|t/'P) at the marked point of the power curve. 


These results hold not only for the cubic nonlinearity, but also for a much wider 
class of nonlinearities. While it has been known that T^T-symmetric dissipative 
systems share some important properties with conservative systems, the results 
in this paper reveal that certain types of non-PT-symmetric dissipative systems 
can also share such properties of conservative systems (such as the existence of 
soliton families). 

Our results also shed light on a more general question: what non-V'T- 
symmetric complex potentials in the NLS-type equations o and (I93|) admit 
continuous families of solitons? In the absence of T^Tsymmetry, the existence 
of a constant of motion in the stationary soliton equation is critical for the ex¬ 
istence of soliton families. We have shown that such a constant of motion exists 
only for special potentials of the form ([5|) , assuming this constant of motion is a 
continuous deformation of that from the potential-free equation. Since this as¬ 
sumption is reasonable, we conjecture that the only non-VT-symmetric complex 
potentials which admit soliton families are those of the special form W- 

It should be pointed out that the question of solitons in non-PT-symmetric 
potentials is closely related to the question of non-PT-symmetric solitons in 
T^T-symmetric potentials. Indeed, for T^T-symmetric potentials of the same 
special form ®, where g{x) is taken to be even, it has been shown numerically 
that symmetry breaking of solitons can occur m- As a consequence, continuous 
families of non-'PT-symmetric solitons exist in a VT -symmetric potential. This 
symmetry breaking is surprising since it is forbidden in generic PT-symmetric 
potentials [28] . Analytical understanding of this symmetry breaking is still an 
open question, however, based on the analysis in this paper, it is hopeful that 
this symmetry breaking can now be analytically studied. But this lies outside 
the scope of the present article. 

In the end, we mention that bifurcation of soliton families from linear modes 
occurs in special forms of two-dimensional non-PP-symmetric complex poten¬ 
tials as well [29]. Analytical understanding of such bifurcations in two spatial 
dimensions is a more challenging question which merits further investigation. 
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